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Abstract: We study the critical behaviour of the 2d dodecahedron spin model and 
investigate the conjecture that the discrete model describes the same continuum theory 
as the 0(3) non-linear sigma model. In particular, we found that the anisotropy of 
the magnetization A(z) measured in a fixed physical volume decreases with increasing 
correlation length, at least up to £ ~ 1000. 


One of the most studied quantum field theories is the 2d non-linear sigma model. The 
lattice regularized version for the 0(3) case is given by the action 

A(S) =--^2S(x)S(x + fi) (1) 

where S(x) is a unit 0(3) vector, S 2 (x) = 1. The model is widely believed to be asymp¬ 
totically free (AF) and to have a dynamically generated mass gap. 

Here we discuss the analogous model when the spin variables are allowed to take only 
discrete values pointing towards the vertices of a Platonic regular solid, e.g. icosahedron 
(I) or dodecahedron (D). Accordingly, the symmetry of the discrete model is reduced to a 
discrete subgroup V of 0(3) (in the actual case, to the icosahedral group with 60 elements). 
At small g (low temperature) the discrete symmetry is spontaneously broken while at large 
g the system is in the symmetric phase under V. Numerical evidence shows that there is 
a second order phase transition at some g = g c > 0. 

A few years ago Patrascioiu and Seiler jlj] observed that within their statistical errors 
the physical results for the dodecahedron model for g \ g c are consistent with those of the 
0(3) model and conjectured that the two models are equivalent in the continuum limit. 
The authors saw this observation as an argument supporting their unorthodox view that 
the 0(3) model is not AF. 

The numerical evidence of [||] inspired us to study this question. In j2] we pointed 
out that the equivalence is not in contradiction with the 0(3) model being AF. Of course, 
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Figure 1: Deviation of gn(z) from the 0(3) Figure 2: The FSS function £(2 L)/£(L) vs 
result vs. 1/£(L). 1/z = £(L)/L. 


for the discrete model this is not a statement about the bare coupling (which goes to 
a non-zero constant in the continuum limit) but about some physical running coupling, 
say the Liischer-Weisz-Wolff coupling g(L) [|3j]. In pj] we measured the renormalized zero- 
momentum 4-point coupling gn and have found that the I and D models approach the 
0(3) model with increasing £. In that work we went up to correlation length 1 * £ ~ 300 and 
reached precision of 0(0.1%). 

In H| we also compared the finite size scaling (FSS) function of the LWW coupling in 
the discrete models to that of the 0(3) model, with a similar conclusion. 

Recently Caracciolo et al. 0 argued that the operator which has to be added to 
the 0(3) action to break the original 0(3) symmetry down to V is a relevant operator and 
concluded that the equivalence cannot hold. They suggested that the numerical evidence is 
misleading and the discrete model should belong to a different universality class, although 
they expected this to show up only at £ > 200. At this conference they have presented 
measurements of the FSS function and concluded that the discrete model starts to depart 
from the 0(3), although at much larger correlation length, £ ~ 10 5 [||. We shall reflect on 
their arguments in the concluding section. 

For the discrete models we have measured gn(z), the renormalized coupling at a fixed 
finite physical volume, z = L/£(L), where £(L) is the “second-moment correlation length”. 
Fig- m shows the deviation of gR,(z) from 0(3) at z = 2.32 as a function of 1/£(L), both for 
the I and D case j2[]. 

The FSS function £(2L)/£(L) for the icosahedron and dodecahedron models [|j is 
shown in fig.|^. Larger symbols refer to larger correlation length, the largest one being 
£ ss 300. The 0(3) curve is taken from Q. The comparison shows that the icosahedron 
and dodecahedron models approach the 0(3) result with increasing correlation length. 
(The results for tetrahedron, octahedron and cube - not shown here - depart from it.) 

The restoration of the 0(3) symmetry can be tested, however, in a more direct way by 
measuring the direction of average magnetization in a fixed physical volume and checking 
whether its distribution becomes uniform in the continuum limit or still prefers the original 

1 Here and below by £ we denote the infinite volume correlation length. When not measured directly it 

is calculated from the FSS function of 
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discrete directions. To do this we introduce the quantity 


In( S) 



( 2 ) 


where S is an arbitrary 0(3) vector, R runs over the set of vectors of the actual discrete 
model, {Rj, i = 1,..., N} while R' runs over the analogous set of the dual regular solid 
(i.e. the unit vectors pointing towards the face centers of the original regular solid). For 
example, in the case of the dodecahedron N = 20 and the dual model is the icosahedron 
with N' = 12. The value of n is chosen to be the smallest integer for which J n (S) is not 
identically zero. For the dodecahedron this is n = 6. Obviously, the integral of I n ( S) over 
S is zero, and I n { S) preserves the discrete symmetry V. We fix the normalization by the 
condition I n ( R) = 1 (where R is one of the discrete directions). Up to the normalization 
constant our I n ( S) coincides with the one introduced by Caracciolo et al. Q. 2 

Consider now a fixed physical volume given by z = L/£(L) and define the anisotropy 

A (z) = » wh ere M = ^ R(z) . (3) 


This quantity is a direct measure of the possible restoration of the 0(3) symmetry in the 
continuum limit. Fig. 3 shows the anisotropy A(z) for the dodecahedron at z = 1.888, 1.668 
and 1.525 against l/£. As the plot shows, A(z) decreases with increasing £, at least up to 
our last value, £ ~ 1000. This means that the dodecahedron model approaches the 0(3) 
model. If the decrease persists and lim^ 00 A(z;£) = 0 then they are equivalent in the 
continuum limit. 

Similarly to Caracciolo et al. Ji|] we also considered the mixed model 

A(S) = --VS(,)S(, + A) - h V/ n (S(x)), (4) 

q z — J z —' 

X,fl X 


which interpolates between the 0(3) model (h = 0) and the discrete model ( h = oo). 
Here one expects a critical line g c {h) connecting <? c (0) = 0 with the critical point g c {oo) 
of the discrete model. For g < g c (h) the discrete symmetry T> is spontaneously broken, 
while for g > g c {h ) the system is in its symmetric phase. Since (by the standard wisdom) 
the correlation length at g > 0, h = 0 is finite, one expects that for h <C m 2 (g ) (where 
m(g) = 1/£(<?) is the inverse of the 0(3) correlation length) no phase transition can occur. 
This is because the effect of the corresponding term on the independently fluctuating 
regions of size £(<?) is negligible. (The situation here is in sharp contrast to the 0(2) case 
where in the massless phase an arbitrarily small external field can have a drastic effect on 
the system in an infinite volume.) 

One can study the RG flow in the ( g , h) plane. The scenario advocated by Caracciolo 
et al. Q means that the RG flow along the critical line g c (h ) is upwards, towards the h = oo 
point, while if the two models are in the same universality class the flow should be directed 

2 Although this is not directly seen since they used a concrete representation for the discrete spin vectors. 


- 3 - 


PrHEP hep2001 





International Europhysics Conference on HEP 



Figure 3: The anisotropy A(z) vs. £. Figure 4: RG flow in the (g, h ) plane. On 

the vertical axis the values of h/{h + 1) are 
plotted. 


downwards, towards g = h = 0. A decisive answer would be provided by measuring the 
flow in the vicinity of the critical line. For example, consider a point (g, h ) corresponding to 
a fixed physical volume, z = L/£(L) and with anisotropy A(z). To determine the direction 
of the flow one has to find a point ( g',h') and size L' such that f'{L') = s£(L), L' = sL 
(i.e. z' = z) and A'(z) = A(z), with a scale factor s < 1. At present we have data only for 
a few points, with s = 1/2, at z ~ 1.9. The matching pairs (g, h ) —> (g 1 , h') are shown in 
fig. ^ where the values of L and L' are also indicated. Accordingly, our largest correlation 
length was £(L) ~ 60 which corresponds to £ ~ 130. (The dashed line representing g c (h) is 
drawn in an ad hoc way, the value /3 c (oo) = l/g c (oo) ~ 2.15 is taken from j7j.) The pattern 
in fig. H suggests a downwards flow, but one would need points at larger correlation length 
to get a convincing answer. 

Now we would like to discuss the scenario suggested by Caracciolo et al.||, I). They 
expand the correlation functions (defined in a flnite L x L volume) in double power series 
in g and h. Although the coefficients of this expansion diverge as L —> oo the functions 
appearing in the corresponding RG equation stay flnite in this limit. The assumption they 
implicitly make is that the RG equation obtained remains valid if both g and h are small 
enough independently. Below we shall argue that the RG flow obtained by PT is valid 
only in a small wedge around the h = 0 axis, for h <C 1 / R 2 where R = min{£(< 7 ), L}. In 
the extreme case of g = 0, L = finite (and for simplicity with the usual magnetic field 
coupled to S z (x )) one has (S z ) = 1/tanh (hL 2 ) — l/(hL 2 ). For fixed L this has a power 
series expansion in h with infrared divergent coefficients. This expansion becomes, however, 
invalid for hL 2 1 where the true answer is (S z ) ~ 1 — 1 /(hL 2 ). A similar behaviour 
is expected for finite g, with L replaced by min{£(g),L}. The reason for the existence of 
these two regimes is obvious: the system behaves differently depending on whether the 
external field is strong enough to keep the (independently fluctuating) regions of size R in 
the preferred direction or not. 

The argument given above is also valid for an hl n ( S) term in the action. There is, 
however, an important difference: in the latter case the remnant discrete symmetry is 
spontaneously broken at h > h c (g). This breaking is associated with tunneling between 
different minima of —hl n { S), i.e. the dodecahedron directions. It is quite obvious that this 
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tunneling cannot be described by PT: one can replace —I n ( S) by a function which has one 
single minimum in one of the preferred discrete directions and has the same perturbative 
expansion around this minimum as — J n (S). PT will not notice any difference between the 
two cases while in the second case no discrete symmetry is left hence no phase transition 
will occur by increasing h. Although these are only speculative arguments, we think they 
could be made more rigorous. 

As mentioned above, in |5[] the authors also found that the FSS function for the do¬ 
decahedron starts to show a deviation from the 0(3) model at huge correlation length, 
£ ~ 10 3 4 5 6 7 8 . Our criticism is the same as given by Patrascioiu and Seiler ||]: 1) for £ ~ 10 5 
and L = 150 the physical size of the system is practically zero and this small distances are 
dominated still by cutoff effects, 2) at such huge estimated £ one cannot be sure that the 
system is in the symmetric phase. 

Summarizing, we think that the spontaneous breaking of the discrete symmetry is 
a non-perturbative critical phenomenon, even in the vicinity of the g = h = 0 point. 
Concerning the numerical evidence - in particular the anisotropy measurements up to 
£ ~ 1000 and the fact that in this range of £ Caracciolo et al. do not see either a deviation 
in the FSS behaviour - indicate that the dodecahedron (and the icosahedron) models 
describe the same continuum theory as the 0(3) model. The question, however, deserves 
further numerical study. 

The authors thank Janos Balog, Sergio Caracciolo, Martin Liischer, Andrea Montanari, 
Andrea Pelissetto and Peter Weisz for useful discussions and correspondence. This work 
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Human Potential Programme under HPRN-CT-2000-00145 Hadrons/Lattice QCD, BBW 
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